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Abstract 

After reexamining the above barrier diffusion problem where we notice that the wave packet col- 
hsion implies the existence of multiple reflected and transmitted wave packets, we analyze the way 
of obtaining phase times for tunnehng/reflecting particles in a particular colliding configuration 
where the idea of multiple peak decomposition is recovered. To partially overcome the analytical 
incongruities which frequently rise up when the stationary phase method is adopted for computing 
the (tunneling) phase time expressions, we present a theoretical exercise involving a symmetrical 
collision between two identical wave packets and a unidimensional squared potential barrier where 
the scattered wave packets can be recomposed by summing the amplitudes of simultaneously re- 
flected and transmitted wave components so that the conditions for applying the stationary phase 
principle are totally recovered. Lessons concerning the use of the stationary phase method are 
drawn. 
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The recent developments of nanotechnology brought new urgency to study the tunnehng 



time as it is directly related to the maximum attaina 
In parallel, a series of recent experimental results 



e speed of nanoscale electronic devices. 



flB.Q 



some of them corroborating 



with the possibility of superluminal tunneling speeds for photons, have revived an interest 
in the tunnelling time analysis B,flQB,a. on the t™.o„t, people Mvetned to 
introduce quantities that have the dimension of time and can somehow be associated with 
the passage of the particle through the barrier or, strictly speaking, with the definition of the 
tunneling time. Since a long time these efforts have led to the introduction of several time 



definitions 



ime. bmce a long time tnese en 

0, 0, Q, Q, HQ, Q, H, 0: 



191], some of which are completely unrelated 



to the others, which can be organized into three groups: (1) The first group comprises a 
time-dependent description in terms of wave packets where some features of an incident 
wave packet and the comparable features of the transmitted packet are utilized to describe 
a delay as tunneling time Q]. (2) In the second group the tunneling times are computed 
with basis on averages over a set of kinematical paths, whose distribution is supposed to 
describe the particle motion inside a barrier, i. e. Feynman paths are used like real paths 
to calculate an average tunneling time with the weighting function exp [i S x{t)/h], where 
S is the action associated with the path x{t) (where x(t) represents the Feynman paths 
initiated from a point on the left of the barrier and ending at another point on the right of 
it 2^). The Wigner distribution paths [0], and the Bohm approach 21, 22] are included 
in this group. (3) In the third group we notice the introduction of a new degree of freedom, 
constituting a physical clock for the measurements of tunnelling times. Separately, it stands 
by itself the dwell time approach. The time related to the tunneling process is defined by 
the interval during which the incident flux has to exist and act, to provide the expected 

Q n 

accumulated particle storage, inside the barrier [7[. The methods with a Larmor clock [11[ 
or an oscillating barrier 23|] are comprised by this group. 



There is no general agreement [sl, Q] among the above definitions about the meaning of 
tunneling times (some of the proposed tunneling times are actually traversal times, while 
others seem to represent in reality only the spread of their distributions) and about which, if 
any, of them is the proper tunneling time, in particular, due to the following reasons [5]: (a) 
the problem of defining tunnelling times is closely connected with the more general definition 
of the quantum-collision duration, and therefore with the fundamental fact that in quantum 
mechanics, time enters as a parameter rather than an observable to which an operator can 
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be assigned (b) the motion of particles inside a potential barrier is a quantum phenomenon, 
that till now has been devoid of any direct classical limit; (c) there are essential differences 
among the initial, boundary and external conditions assumed within the various definitions 
proposed in the literature; those differences have not been sufficiently analyzed yet. In 
particular, the study of tunneling mechanisms is embedded by theoretical constructions 
involving analytically-continuous gaussian, or infinite-bandwidth step pulses to examine the 
tunneling process. Nevertheless, such holomorphic functions do not have a well-defined 
front in a manner that the interpretation of the wave packet speed of propagation becomes 
ambiguous. Moreover, infinite bandwidth signals cannot propagate through any real physical 
medium (whose transfer function is therefore finite) without pulse distortion, which also leads 
to ambiguities in determining the propagation velocity during the tunneling process. For 
instance, some of the barrier traversal time definitions lead, in tunneling time conditions, to 
very short times which can even become negative, precipitately inducing to an interpretation 
of violation of simple concepts of causality. Otherwise, negative speeds do not seem to create 
problems with causality, since they were predicted both within special relativity and within 
quantum mechanics 171]. A possible explanation of the time advancements related to the 
negative speeds can come, in any case, from consideration of the very rapid spreading of the 
initial and transmitted wave packets for large momentum distribution widths. Due to the 
similarities between tunnelling (quantum) packets and evanescent (classical) waves, exactly 
the same phenomena are to be expected in the case of classical barriers (we can mention 
the analogy between the stationary Helmholtz equation for an electromagnetic wave packet 
- in a waveguide, for instance - in the presence of a classical barrier and the stationary 
Schroedinger equation, in the presence of a potential barrier , |2J] ) . The existence of 



such negative times is predicted by relativity itself, on the basis of the ordinary postulates 
jsl, and they appear to have been experimentally detected in many works 
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In this extensive 



phase times [27 



. y explored scenario, the first group quoted above contains the so-called 
which are obtained when the stationary phase method (SPM) {29! is 
employed for obtaining the times related to the motion of the wave packet spatial centroid 

which adopts averages over fluxes pointing in a well-defined direction only, and has recourse 

j 

to a quantum operator for time [5| . Generically speaking, the SPM essentially ena 
parameterize some subtleties of several quantum phenomena, such as tunneling 



es us to 

,y,[i6i, 

resonances 30|], incidence-reflection and incidence-transmission interferences [311] as well as 
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the Hartman Effect 



34| and its superluminal traversal time interpretation |5|,l7|,ll8|. In fact, 



Bfl 



it is the simplest and commonest approximation method for describing the group velocity of 
a wave packet in a quantum scattering process represented by a collision of a particle with 



a potential barrier 



28 



34 



35]. 



Our attention is particularly concentrated on some limitations on the use of the SPM for 
deriving tunneling times for which we furnish an accurate quantification of the analytical 
incongruities which restrict the applicability of this method. We introduce a theoretical 
construction involving a symmetrical collision with a unidimensional square potential where 
the scattered wave packets can be reconstructed by summing the amplitudes of the reflected 
and transmitted waves in the scope of what we denominate a multiple peak decomposition 
analysis js^ in a manner that the analytical conditions for the SPM applicability are totally 
recovered. 

Generically speaking, the SPM can be successfully utilized for describing the movement 
of the center of a wave packet constructed in terms of a momentum distribution g{k — k^) 
which has a pronounced peak around k^. By assuming the phase which characterizes the 
propagation varies sufficiently smoothly around the maximum of g{k — ko), the stationary 
phase condition enable us to calculate the position of the peak of the wave packet (highest 
probability region to find the propagating particle). With regard to the tunneling effect, the 
method is indiscriminately applied to find the position of a wave packet which traverse a 
potential barrier. For the case we consider the potential barrier 



V{x) 



Vn 



X G [-L/2, L/2] 



(1) 



X ^ [-L/2, L/2] 

it is well known that the transmitted wave packet solution (x > L/2) calculated by means of the 



Schroedinger formalism is given by 

dk 



36] 



2tt 



g{k - ko) \T{k,L)\ exp 



ik(x- L/2) -i t + i e(k, L) 

2 m 



where, in case of tunneling, the transmitted amplitude is written as 

„4 ^ -I 



\nk,L)\ 



1 + 



w 



Ak^ p^{k) 

and the phase shift is obtained in terms of 



sinh^ [p{k) L] 



e{k,L) 



arctan 



2 k' 



w 



kp{k) 



tanh [p{k) L] } , 



(2) 



(3) 



(4) 



for which we have made exphcit the dependence on the barrier length L and we have adopted 
p{k) = [w"^ — k"^)^ with w = (2mV"o)^ and h = 1. Without thinking over an eventual 
distortion that |T(A;,L)| causes to the supposedly symmetric function g{k — k^), when the 
stationary phase condition is applied to the phase of Eq. ([3]), we obtain 



±Sk{x-L/2)-^t + e{k,L) 
ak I 2 m 







0. (5) 



m dk 

The above result is frequently adopted for calculating the transit time of a transmitted 
wave packet when its peak emerges at x = L/2, 

m de{k,a^j^^) 



^max dk 

2m L sinh (a) cosh (a) - (2 kj^^ - w^) fc^^ a \ 

a \ 4 kl^^ {w' - kiJ + sinh^ (a) J ^ ^ 

where we have introduced the parameter a = [w"^ — A;^^^)^ L. The concept of opaque limit 



appears when we assume that fc^ax is independent of L and then we make a tends to oo |l8 |. 

In this case, the transit time can be rewritten as 

OL 2 m 

= k oik V 

In the literature, the value of fc^ax is frequently approximated by k^, the maximum of g{k—ko), 
which, in fact, does not depend on L and, under the theoretical point of view, could lead to 
the superluminal transmission time interpretation ,|4Qi. 



It would be perfectly acceptable to consider k^^^ = k^ for the application of the stationary 
phase condition if the momentum distribution g{k — ko) centered at kg had not been modified 
by any boundary condition. This is the case of the incident wave packet before colliding 
with the potential barrier. Our criticism concern with the way of obtaining all the above 
results for the the transmitted wave packet. It has not taken into account the bounds and 
enhancements imposed by the analytical form of the transmission coefficient. 

To perform the correct analysis, we should calculate the right value of k^^^ to be substi- 
tuted in Eq. (jH]) before taking the opaque limit. We should be obliged to consider the relevant 
amplitude for the transmitted wave as the product of a symmetric momentum distribution 
g{k — ka) which describes the incoming wave packet by the modulus of the transmission am- 
plitude T[k, L) which is a crescent function of k. The maximum of this product representing 



the transmission modulating function would be given by the solution of the equation 



g{k-K) \T{k,L) 



g'{k-K) ^\T{k,L) 



0. (8) 



,gik-h) \T{k,L) 

Obviously, the peak of the modified momentum distribution is shifted to the right of k^ so 
that k^^^ has to be found in the interval ]fco, w[. It could be numerically confirmed that fc^ax 
presents an implicit dependence on L so that, by increasing the value of L with respect to 
a, the value of k^^^ to be utilized in Eq. ([6]) would increase until L reaches certain values 
for which the modified momentum distribution becomes unavoidably distorted. In this case, 
the relevant values of k are concentrated around the upper boundary value w. 

Due to the filter effect, the amplitude of the transmitted wave is essentially composed 
by the plane wave components of the front tail of the incoming wave packet which reaches 
the first barrier interface before the peak arrival. Meanwhile, only whether we had cut the 
momentum distribution Oj^at a value of k smaller than w, i. e. ~ {l—6)w, the superluminal 
interpretation of the transition time could be recovered. In this case, independently of 
the way as a tends to oo, the value assumed by the transit time would be approximated by 
t^ ^ 2m/w6 which is a finite quantity. Such a finite value would confirm the hypothesis 
of superluminality. However, the cut off of the momentum distribution at k ^ {1 — 5)w 
increases the amplitude of the tail of the incident wave and it contributes so relevantly as 
the peak of the incident wave to the final composition of the transmitted wave so that an 
ambiguity in the definition of the arrival time is created. 

We are particularly convinced that the use of a step-discontinuity to analyze signal trans- 
mission in tunneling processes deserves a more careful analysis than the immediate applica- 
tion of the stationary phase method since we cannot find an analytic-continuation between 
the above barrier case solutions and the helow barrier case solutions. A suggestive possibility 
can ask for the use of the multiple peak decomposition technique developed for the above 
barrier diffusion problem 35|]. Thus, in a similar framework, we suggest a suitable way for 
comprehending the conservation of probabilities for a very particular tunneling configuration 
where the asymmetric aspects discussed up to now can be totally eliminated, with the phase 
times being accurately calculated. By means of such an experimentally verifiable exercise, 
we shall be able to understand how the filter effect can analytically affect the calculations 
of transit times in the tunneling process. 

In order to recover the scattered momentum distribution symmetry conditions for cor- 
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rectly applying the SPM, we assume a symmetrical colliding configuration of two wave 
packets traveling in opposite directions. By considering the same barrier represented in (Q, 
we solve the Schroedinger equation for a plane wave component of momentum k for two iden- 
tical wave packets symmetrically separated from the origin a; = 0. At time t = —{rriL) / {2ko) 
chosen for mathematical convenience, we assume they perform a totally symmetric simulta- 
neous collision with the potential barrier. The wave packet reaching the left(right)-side of 
the barrier is thus represented by 

^IjMr) (^x, t)= dk g{k - ko)(l)^^''^ {k, x) exp [-i E t] (9) 

where, as a first approximation, we are assuming that the integral can be naturally extended 
from the interval [0, w] to the interval [0, oo\. Its range of validity can be controlled by the 
choice of the width A/c of the momentum distribution g{k — ko) (with ko > 0), i. e. Ak has to 
be bounded by the barrier high (V^) in order to avoid the contribution of the above-barrier 
frequencies (or energies) contained in the considered wave packet (which eventually become 
important as the tunnelling components are progressively damped down). 

By assuming that ^^^''^k, x) are Schroedinger equation solutions, when the wave packet 
peaks simultaneously reach the barrier ( at the time t = — (mL)/(2A;„)) we can write 

ih^''\k,x) = exp[±ikx] + R^^''\k,L)e^p[Tikx] x < -L/2{x > L/2), 

^5^(^)(A;,x) = a^^''\k)exp[^px]+ (3l^''\k)exp[±px] - L/2 < x < L/2, 
pY''\k,x) = Tl'^''\k,L)exp[±ikx] X > L/2{x < -L/2). 

where the upper(lower) sign is related to the index L{R). By assuming the conditions for the 

continuity of 0^'^ and their derivatives at x = —L/2 and x = L/2, after some mathematical 

manipulations, we can easily obtain 

pL,B,, r .,^,J exp[i9(fc,L)] [l-exp[2p(fc)L]] ) 

{k, L) = exp [-Z ^ L] I ^ _ ^^^^ ^^^^^^ I (10) 

and 

where Q{k,L) is given by the Eq. (jll) and I^{k,L) and T^{k,L) as well as R%{k,L) and 
T]^{k,L) are intersecting each other. By analogy with the procedure summing amplitudes 
which we have adopted in the multiple peak decomposition scattering [35|, such a pictorial 
configuration obliges us to sum the intersecting amplitude of probabilities before taking their 



0l^W(A;,x) 
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squared modulus in order to obtain 

R\^-{k^L)+T-^(KL) = eM-^kL] ^ ^xp [p(A:) L] + exp e(fc, L)] 



1 + exp [p{k) L] exp [i Q{k, L)] 
exp{~i[kL + ip{k,L)]} (12) 



with 



/, f 2kp(k) smh\p(k)L] ] 

A;, L = arctan <^ — '^^^^ }. (13) 

The important information we get from the relation given by f|T2|) is that, differently from 
the previous standard tunneling analysis, by adding the intersecting amplitudes at each side 
of the barrier, we keep the original momentum distribution undistorted since \R^'^{k, L) + 
Tg'^{k, L) \ is equal to one. At this point we recover the most fundamental condition for the 
applicability of the SPM which allows us to accurately find the position of the peak of the 
reconstructed wave packet composed by reflected and transmitted superposing components. 

The phase time interpretation can be, in this case, correctly quantified in terms of the 
analysis of the new phase tf{k, L). By applying the stationary phase condition to the recom- 
posed wave packets, the maximal point of the scattered amplitudes g{k — ko)\R^'^{k, L) + 
Tg'-^ (k, L)\ are accurately given by k,^^^ = kg so that the traversal/reflection time or, more 



generically, the scattering time, results in 

^ _ m dip{k,a^^) 

tn 



dk 



2m L sinh (q) — a 

(14) 



k(,a 2 kl — w'^ + w'^ cosh (a) 

with a previously defined. It can be said metaphorically that the identical particles repre- 
sented by both incident wave packets spend a time of the order of ti^ inside the barrier before 
retracing their steps or tunneling. In fact, we cannot differentiate the tunneling from the 
reflecting waves for such a scattering configuration. The point is that we have introduced 
a possibility of improving the efficiency of the SPM in calculating reflecting and tunneling 
phase times by studying a process where the conditions for applying the method are totally 
recovered, i. e. we have demonstrated that the transmitted and reflected interfering ampli- 
tudes results in a unimodular function which just modifies the envelop function g{k — k^) by 
an additional phase. The previously appointed incongruities which cause the distortion of 
the momentum distribution g{k — /cq) are completely eliminated in this case. At the same 
time, one could argue about the possibility of extending such a result to the standardly 
established tunneling process. We should assume that in the region inside the potential 
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barrier, the reflecting and transmitting amplitudes should be summed before we compute 
the phase changes. Obviously, it would result in the same phase time expression represented 
by ([II]). In this case, the assumption of there (not) existing interference between the mo- 
mentum amplitudes of the reflected and transmitted at the discontinuity points x = —L/2 
and X = L/2 is purely arbitrary. Consequently, it is important to reinforce the argument 
that such a possibility of interference leading to different phase time results is strictly re- 
lated to the idea of using (or not) the multiple peak (de) composition in the region where 
the potential barrier is localized. 




FIG. 1: Time ratios for the standard tunneling and the new scattering process. The ratios R{a) 
and R't'ia) can be understood as transmitted times in the unities of the classical propagation time 
r. Both of them present the same asymptotic behavior which recover the theoretical possibility 
of a superluminal transmission in the sense that by now, under the point of view of the analytical 
limitations, the SPM can be accurately applied. 

In order to illustrate the difference between the standard tunneling phase time tx and the 
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the alternative scattering phase time we introduce the new parameter n = kl^^^/ and 
the classical traversal time r = (mL) / k^-^^^ in order to define the ratios 



, , It 2 f cosh (a) sinh (a) — Q n (2n — 1) 
Knp[a) = — = - <^ - 



a 



and 



[4n (1 — n) + sinh^ (a)] 

/ ^ 2 f n a + sinh (a) 

' T a [2n-l + cosh (a) 



(15) 



(16) 



which are plotted in the Fig.([T]) for some discrete values of n varying from to 1, from which 
we can obtain the commonest limits given by 



lim {Rria)} = lim {Rria)} = 



and 



lim {Rria)} = 1 + ^, 
a^O In 



hm {RUa)] 



1 + 



1 



n 



(17) 



(18) 



Both of them present the same asymptotic behavior which, at first glance, recover the 
theoretical possibility of a swj»er/wmma/ transmission in the sense that, by now, the SPM can 
be correctly applied since the analytical limitations are accurately observed. At this point, it 
is convenient to notice that the superluminal phenomena, observed in the experiments with 
tunnelling photons and evanescent electromagnetic waves , has generated a lot 

of discussions on relativistic causality. In fact, superluminal group velocities in connection 
with quantum (and classical) tunnelings were predicted even on the basis of tunneling time 
definitions more general than the simple Wigner's phase-time 28|] (Olkhovsky et ai, for 
instance, discuss a simple way of understanding the problem [s]). In a causal manner, it 
might consist in explaining the superluminal phenomena during tunnelling as simply due to 
a reshaping of the pulse, with attenuation, as already attempted (at the classical limit) 32| . 
i. e. the later parts of an incoming pulse are preferentially attenuated, in such a way that 
the outcoming peak appears shifted towards earlier times even if it is nothing but a portion 
of the incident pulse forward tail 0, 33|. In particular, we do not intend to extend on the 
delicate question whether superluminal group-velocities can sometimes imply superluminal 
signalling, a controversial subject which has been extensively explored in the literature about 
the tunneling effect (js] and references therein). 

Turning back to the scattering time analysis, we can observe an analogy between our 
results and the results interpreted from the Hartman Effect (HE) analysis 3J]. The HE 
is related to the fact that for opaque potential barriers the mean tunnelling time does not 
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depend on the barrier width, so that for large barriers the effective tunnelhng-velocity can 
become arbitrarily large, where it was found that the tunnelling phase-time was independent 
of the barrier width. It seems that the penetration time, needed to cross a portion of a bar- 
rier, in the case of a very long barrier starts to increase again after the plateau corresponding 
to infinite speed proportionally to the distance 4l|. Our phase time dependence on the 
barrier width is similar to that which leads to Hartman interpretation as we can infer from 
Eqs. fl 17111 81) . Only when a tends to we have an explicit linear time-dependence on L given 
by 

2 m L / 1 \ , ^ 

■w \ n J 

which agree with calculations based on the simple phase-time analysis where tx = 
(l + 2^). However, it is important to emphasize that the wave packets for which we 
compute the phase times illustrated in the Fig.([T]) are not effectively constructed with the 
same momentum distributions. The phase <d{k,L) appears when we treat separately the 
momentum amplitudes g{k — k^) \T{k, L)\ and g{k — ko)\R{k, L)\ and the other one ip{k, L) 
appears when we sum the amplitudes g{k — ko) \T{k, L) + R{k, L)\ = g{k — k^) in order to 
obtain a symmetrical distribution which thus "requalify" the SPM for exactly determine the 
time-position of the peak of a wave packet. In this sense, as a suggestive possibility for par- 
tially overcoming the incongruities (here appointed and quantified) which appear when we 
adopt the SPM framework for obtaining tunneling phase times, we have claimed for the use 
of the multiple peak decornposition [35| technique presented in the study of the the above 
barrier diffusion problem 35|]. We have essentially suggested a suitable way for compre- 
hending the conservation of probabilities for a very particular tunneling configuration where 
the asymmetry presented in the previous case was eliminated, and the phase times could 
be accurately calculated. An example for which, we believe, we have provided a simple but 
convincing resolution. 

In a more extended context, there also have been some attempts of yielding complex 
time delays, ultimately due to a complex propagation constant. This has caused some 
controversies with denying the physical reality to an imaginary time 3^. In parallel to 
the most sensible candidate for tunneling times , a phase-space approach have been 
used to determine a semi-classical traversal time 37|. This semi-classical method makes 
use of the concept of complex trajectories which, by its turn, enables the definition of 
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real traversal times in the complexified phase space. It is also commonly quoted in the 
context of testing different theories for temporal quantities such as arrival, dwell and delay 
times and the asymptotic behavior at long times [18|, |38(]. The configuration we have 
introduced in this manuscript suggests that perhaps the idea of complexifying time should 
be investigated for some other scattering configurations. We let for a subsequent analysis 
the suggestive possibility of investigating the validity of our approach when confronting with 
the phenomenon o 
potential barriers 



one- dimensional non-resonant tunnelling through two successive opaque 



391] and with the intriguing case of multiple opaque barriers [40|]. Still 
concerning with the future theoretical perspectives, the symmetrical colliding configuration 
also offers the possibility of exploring some applications involving soliton structures. Finally, 
we believe all the above arguments reinforce the assertion that it is necessary to continue 
the search for a generalized framework where barrier traversal times can be computed. 
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